Impact of Auger processes on carrier dynamics in graphene 
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The linear, gapless bandstmcture of graphene provides ideal conditions for Auger processes. They are of 
great importance for fundamental research and technological applications, since they qualitatively change the 
canier dynamics. Time-, momentum-, and angle-resolved microscopic calculations reveal an efficient impact 
excitation giving rise to a significant multiplication of optically excited carriers and a remarkable Coulomb- 
induced carrier cooling effect. We present an analytic expression for the carrier multiplication predicting ideal 
conditions for its still missing experimental observation. 



Graphene consists of a single layer of carbon atoms ar- 
ranged in a honeycomb crystal lattice. Its reduced dimen- 
sionality and the unusual linear, gapless bandstmcture leads 
to unique electronic and optical properties suggesting its ap- 
plication in various opto-electronic devices yj-yl]. The key 
prerequisite for the realization of such applications is a better 
microscopic understanding of the teniporal evolution of the 
carrier density in a non-equilibrium |j2|]. In this context, it 
is of great importance to study the efficiency of Auger scatter- 
ing channels, which change the carrier density by bridging the 
valence and the conduction band. In contrast to semiconduc- 
tors having a bandgap and a parabolic bandstmcture, graphene 
shows efficient Auger- type scattering 10-01 ■ 

In the last years, a number of theoretical |0-l8l and experi- 
mental studies |@-{l3l were performed with the aim to obtain 
a detailed understanding of the carrier relaxation dynamics in 
graphene. However, only a few of these studies consider the 
role of Auger scattering channels (0, [sl [Tsl. In this work, 
we obtain new insights into the influence of these channels 
on the carrier density, the carrier temperature, and the chem- 
ical potential in optically excited graphene. Based on a mi- 
croscopic approach offering access to time-, momentum-, and 
angle-resolved carrier relaxation dynamics, we predict (i) a 
significant carrier multiplication (CM) even in the presence 
of directly competing phonon-induced recombination, and (ii) 
a counter-intuitive Coulomb-induced carrier cooling effect, a 
process widely overseen in current literature. We further ad- 
dress the efficiency of Auger-type channels by deriving an 
analytic expression for the strength of the CM reflecting its 
dependence on the pump fluence, temperature, and excitation 
energy. This allows us to give the optimal conditions for its 
still missing experimental observation in graphene. 

Carrier multiplication is defined as the ratio between the 
total number of generated and the initially optically excited 
carrier-hole pairs, reflecting the enhancement of the quan- 
tum yield due to many-particle effects. Under certain con- 
ditions, the absorption of a single photon can generate mul- 
tiple electron-hole pairs, cp. Fig. [T] It has already been 
observed experimentally and studied theoretically in quantum 
dots 1 15 . T3l and carbon nanotubes 1 17 - 23] . In graphene how- 
ever, an efficient CM has been predicted theoretically but 
has not yet been directly measured in experiments. Recent 
studies on CVD-grown multilayer graphene already indicate 
that Coulomb-scattering leads to a change in the carrier den- 
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FIG. 1. (a) Illustration of the impact excitation (IE) increasing the 
number of charge carriers. Two electron-hole pairs are generated by 
the absorption of a single photon - a process called carrier multiplica- 
tion (CM), (b) Schematic carrier distribution showing the efficiency 
of IE at the band edge and displaying the possible scattering partners 
provided by the optically excited non-equilibrium distribution. 



sity 1131. 

The process of CM is driven by Auger scattering consisting 
of two inverse processes: impact excitation (IE) and Auger re- 
combination (AR). Figure[T^ illustrates IE, where an optically 
excited electron scatters down to an energetically lower state. 
At the same time, another electron is excited from the valence 
into the conduction band. As a result, we obtain two electron- 
hole pairs by the absorption of a single photon. In contrast, 
the inverse process of Auger recombination gives rise to a de- 
crease in the charge carrier density (reverse arrows in Fig. [TJ)). 

Our time-resolved calculations are based on the graphene 
Bloch-equations |8], allowing for a consistent consideration 
of carrier-carrier and carrier-phonon scattering on a micro- 



scopic footing 112 IL 12211 . The relaxation dynamics is described 
by a coupled set of differential equations for (i) the occupation 
probability in the state k for electrons and holes (A = c, h), 
(ii) the microscopic polarization pk, which is a measure for the 
transition probability between the bands, and (iii) the phonon 
occupation for different optical and acoustic modes j: 
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Besides the transition frequency Awk, the phenomenological 
phonon decay rate 12311 . and the Bose-Einstein distiibu- 
tion rts as the initial condition for the phonon occupation, all 
terms in in Eq. ([TJ are explicitly time-dependent. The Rabi 
frequency Vl"^' describes the coupling to the exciting vector 
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potential, il,^'^ is the corresponding intraband contribution. 
The many-particle interactions are treated within the second- 
order Bom-Markov approximation resulting in emission F^™ 
and absorption rates F^^ for the phonons. For the micro- 
scopic polarization Pk(i) we obtain diagonal 7k (i) and off- 
diagonal dephasing Uu{t) and for the carrier occupation prob- 
ability a Boltzmann-like scattering equation (sl]: 
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In the case of the Coulomb interaction, the probability to scat- 
ter from the state(Ai , k') to (A, k) is given by: 



Ai^A 
k'->k 



E 

A2,A3,q 



V, 



k,k',q 



(1 



A3 \ An 

Pk'+q)Pk+q' 



(3) 



where V, 
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contains the screened Coulomb matrix ele- 



ment as well as a delta-function ensuring energy conservation 
Jit]. The simultaneous conservation of momentum and energy 
restricts Auger-type processes (IE, AR) to parallel scattering 
along the Dirac cone |4], cp. Fig. [T^. The efficiency of im- 
pact excitation (IE) V£"kf q can be understood in terms of two 
sub-processes denoted in Eqs. (|2]i and (|3]l. The primary sub- 
process (Eq. |2]i k' ^ k bridges both bands and prefers an 
abrupt gradient in the occupation, which is given by a ther- 
mal distribution around the Dirac-point , cp. Fig. [TJ). Here, 
the inverse process of Auger recombination (AR) is strongly 
suppressed by Pauli-blocking. The sum in Eq. [3] includes all 
secondary sub-processes k + q — !> k' + q matching in energy 
and momentum. Note that a non-equilibrium carrier distribu- 
tion is needed to initiate Auger-type scattering, cp. also Fig. 
[TJi. As a result, the efficiency of IE and AR is given by the 
corresponding Coulomb matrix elements and the occupation 
probability of all involved states. 

The numerical evaluation of the Bloch equations enables us 
to track the excited carriers on their way to equilibrium. Our 
calculations reveal that the relaxation dynamics is character- 
ized by a significant carrier multiplication, as shown in Fig. |2] 
Since the efficiency of Coulomb-induced scattering processes 
strongly depends on the number of available charge carriers 
in the system, we study the temporal evolution of CM as a 
function of the optically-induced pump fluence Sopt, cp. Fig. 
|2^. The larger the pump fluence, the more scattering partners 
are available and the faster is the relaxation dynamics. As a 
result, a thermalized distribution is reached already within the 
first tens of fs and the impact of Auger processes is restricted 
to a very short time frame. The direct consequence is a sig- 
nificantly reduced CM for large pump fluences. Moreover, 
the carrier-phonon scattering channels reduce CM for two rea- 
sons: (i) the carrier density decreases via phonon-induced re- 
combination processes and (ii) there is a direct competition 
between Auger and phonon processes, where an excited elec- 
tron can either perform Auger scattering or reduce its energ y 
by emitting a phonon. Figure|2]maps the temporally resolved 
interplay between Coulomb- and phonon-induced channels 
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FIG. 2. (a) Temporal evolution of the Auger-induced carrier mul- 
tiplication (CIVI) for different pump fluences at 300 K. The applied 
optical pulse has an photon energy of 1.5 eV and a duration of 10 fs 
corresponding to recent experiment ||9||. CM is efficient and stable on 
a picosecond timescale at low pump fluences. (b) Time-dependent 
CM at a fixed pump fluence and (c) fluence-dependent of CM at a 
fixed time corresponding to the two straight lines in the contour plot. 



at different pump fluences. Two different regions character- 
ized by CM > 1 (blue area) and CM < 1 (red area) can 
be clearly seen. For a better quantification, two exemplary 
lines of the contour plot are depicted in Fig. |2j3 showing the 
cross-over between both regions. The carrier multiplication 
as a function of time at a fixed pump fluence illustrates the 
efficiency of phonon-induced recombination. On the other 
side, the fluence-dependent cross-over reflects the accelerated 
Coulomb scattering at higher carrier densities. We predict 
a CM of up to a value of 2.5 at low pump fluences. Here, 
CM > 1 remains stable on a ps time scale confirming the 
potential of graphene for applications in opto-electronics. 

Caiiier multiplication can be treated within the quasi- 
equilibrium quantum statistics in consideration of a purely 
Coulomb-induced relaxation dynamics. Energy dissipation 
and recombination due to carrier-phonon scattering reduces 
CM, however, it does not change its qualitative behavior The 
carrier density n and the energy density e of the electronic 
system are determined by the microscopic carrier occupation 
Pk via 
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where and cr^ stand for the spin and valley degeneracy, 
respectively. Furthermore, A is the graphene area and de- 
notes the electronic bandstructure. In equilibrium is given 
by a Fermi distribution and both carrier and energy density 
are functions of the temperature T and the chemical potential 
/i, and vice versa. Therefore, it is possible to completely de- 
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termine the Coulomb-induced final state of the system after 
an optical excitation. In the case of a conventional semicon- 
ductor, where Auger-processes are negligible, the absorbed 
pump fluence Sopt and the optically injected carrier density 
Uopt = '^^optl'fvjJ are both conserved and form new separate 
Fermi distributions in the conduction and the valence band. In 
contrast, in graphene only the energy is conserved, whereas 
the carrier density is a variable quantity giving rise to a sin- 
gle Fermi distribution in both bands. In particular, in the case 
of undoped graphene, the symmetry between electrons and 
holes results in a vanishing chemical potential for the equili- 
brated system. Consequently, in this situation the carrier and 
energy density are determined only by the electronic temper- 
ature. Assuming a linear energy dispersion = ±i^i?|k|, 
Eq. © yields n{T) = ciT^ and e{T) ^ c^T'^ with the co- 
efficients c\ = Trkg/Sup and C2 = 6<^{3)kg/Tri'p, where ( 
denotes the Zeta-function 10] . Starting with an initial thermal 
distribution at the temperature Tq and assuming again an opti- 
cal excitation, we can successively calculate the final state of 
the system / in terms of the energy density £/ = e{To) + eopt, 
the temperature Tf{ef), and the carrier density nf{Tf). Fi- 
nally, we obtain an analytic expression, describing the carrier 
multiplication as a function of the initial temperature Tq, the 
absorbed pump fluency Sopt, and the excitation energy hu: 
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Note that this relation describes the final state of the carrier 
system due to the carrier-carrier scattering and is independent 
of the relaxation time. In particular, Eq. (|5) does not include 
the phonon-induced reduction of CM, which cannot be statis- 
tically treated. Nevertheless, the presented analytic expression 
outlines a systematic estimation of the efficient regime for a 
high carrier multiplication. 

The different dependences of the carrier multiplication 
CM{e opt, To, huj) are illustrated in Figure |3] The compari- 
son between analytic and numerical results shows an excel- 
lent agreement. The highest CM is found for small pump 
fluences, small initial temperatures, and high excitation ener- 
gies. The influence of the pump fluence has already been dis- 
cussed in the previous section - the smaller the pump fluence 
the larger is CM, as shown in Fig. |3^. Now, this dependence 
is quantified by Eq. |5] in particular for small initial temper- 

— 1/3 

atures, CM is proportional to 2^,^/ . For a fixed excitation, 
CM shows only weak dependence of the initial temperature 
To, as displayed in Fig. [Sji- At higher temperatures, it slightly 
decreases resulting from a reduced asymmetry between the 
Auger-type processes of IE and AR. The higher the temper- 
ature, the smoother is the Fermi function and the smaller is 
the suppression of AR due to the Pauli blocking. However, 
compared with the pump fluence, the influence of the initial 
temperature is small in the range between 10 K and 400 K. 
Finally, the appearing CM scales linearly with the excitation 
energy huj, as predicted by the analytic expression and shown 
in Fig. [3};. This can be understood intuitively, since the fur- 
ther away from the Dirac point a carrier is excited, the more 



T 1 1 I I I I I 1 1 1 1 I I I I I 1 1 1 1 I I I I I 



To=300 K,H(0=1.5 eV 
simulation 




b) 




c 


8 r 






1 licat 


6 - 


.& 
■3 


4 - 


S 




U 


2 - 





- 



absorbed pump fluence e^^ [p,Jcm" ] 



^opt-47nJcm" ,HC0=1.5eV 

simulation • 



U— • — 1 * ■ 



_I I I I I I I I I I I I I I I I I I L_ 



c) 



100 200 300 400 

initial temperature Tg [K] 



To=300K, eopt=0.23^J/cm 

• simulation 



■2 




0.5 1 1.5 

excitation energy Hco [eV] 



FIG. 3. Analytically obtained purely Coulomb-induced carrier mul- 
tiplication (red solid line) is compared to the numerical solution from 
Bloch equations (blue dots) for different (a) absorbed pump fluences 
£opt, (b) initial temperatures To, and (c) excitation energies hu). 



IE it can induce on its way back to equilibrium. If the non- 
equilibrium is located close to the Dirac point, i.e. in the case 
of small excitation energies, the scattering probabilities dif- 
fer qualitatively. In contrast to the scattering event depicted 
in Fig. [TJ), both sub-processes can directly be initiated by the 
optically excited carries. This leads to an enhanced efficiency 
for AR, which can even result in CM < 1, cp. Fig. 

In conclusion, ideal condition for an experimental observa- 
tion of CM are a relatively small pump fluence, a low initial 
temperature, and a high excitation energy at the same time. 
The comparison with Fig. |2] shows that the analytically pre- 
dicted efficient regime for CM is also predominant in the pres- 
ence of carrier-phonon scattering. Thus, the presented ana- 
lytic expression gives crucial insights into the dependences of 
CM and helps to estimate its most efficient regime, also un- 
der realistic conditions. Note that deviations from the linear 
bandstructure may reduce the CM. However, we expect non- 
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FIG. 4. Temporal evolution of the (a) temperature and (b) chemical 
potential of an optically excited earner system illustrating the contri- 
bution of Auger-induced processes to carrier cooling and faster equi- 
libration of the chemical potential, respectively. The inset shows the 
momentum-resolved carrier distribution after 100 fs with and with- 
out Auger channels. 

Markov processes to compensate this effect, since they soften 
the strict energy conservation 1I2II1 . 

Another interesting and counterintuitive aspect of the 
Auger processes is a Coulomb-induced carrier cooUng effect. 
In the case of pure carrier-carrier scattering, where the to- 
tal energy of the carrier system is conserved, the electronic 
temperature is a monotonous function of the carrier density. 
Adding carriers to the system via IE, results in a decreased 
temperature. In contrast, AR reduces the number of carriers 
accounting for a carrier heating. We determine the electronic 
temperature as a function of the carrier and energy density, 
which are obtained by the numerical evaluation of the Bloch 
equations including carrier-carrier as well as carrier-phonon 
scattering channels ||24|1 . Figure shows the temporal evo- 
lution of the temperature, illustrating that the predominant IE 
indeed contributes to a more efficient carrier cooling: Switch- 
ing off the Auger channels slows down the decrease in the 
carrier temperature, resulting in a difference of about 400 K 
after 100 fs. The enhanced cooling is also accompanied by a 
faster equilibration of the chemical potential, since the same 
energy is shared by an increasing number of available carriers, 
cp. Fig|4j5. This is further shown by the corresponding carrier 
distributions in the inset. The presence of Auger scattering 
leads to higher occupations at small energies, whereas the ne- 
glect of IE and AR results in a broader carrier distribution. 

In conclusion, we have microscopically investigated the 
role of Auger processes on the ultrafast carrier relaxation dy- 
namics in graphene. We predict a significant carrier multipli- 
cation and an efficient carrier cooling due to the strong impact 
excitation. In spite of the directly competing phonon-induced 
processes, we find a multiplication of carriers by a factor of 
2.5 confirming the potential of graphene as a new material for 
high-efficiency photo-devices. The presented analytic expres- 
sion gives new insights into Auger scattering and guides new 



experiments on graphene and related structures. 
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